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We solve the gap equation for color-superconducting quark matter in the 2SC phase, including 
both the energy and the momentum dependence of the gap, (j) = 4'{ko,'^)- For that purpose a 
complex Ansatz for (j) is made. The calculations are performed within an effective theory for cold 
and dense quark matter. The solution of the complex gap equation is valid to subleading order in 
the strong coupling constant g and in the limit of zero temperature. We find that, for momenta 
sufficiently close to the Fermi surface and for small energies, the dominant contribution to the 
imaginary part of (j> arises from Landau-damped magnetic gluons. Further away from the Fermi 
surface and for larger energies the other gluon sectors have to be included into Imcj). We confirm 
that Im (j) contributes a correction of order g to the prefactor of cf> for on-shell quasiquarks sufficiently 
close to the Fermi surface, whereas further away from the Fermi surface Im^ and Yiecj) are of the 
same order. Finally, we discuss the relevance of Im (j) for the damping of quasiquark excitations. 

PACS numbers: 12.38.Mh, 24.85.-(-p 



I. INTRODUCTION 
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■ Sufficiently cold and dense quark matter is a color superconductor P,|2)S0JS@jI3jI3- In the limit of asymptotically 
^ii^ large quark chemical potentials, ^ Aqcd, quarks are weakly coupled and interact mainly via single-gluon 
QQ . exchange. In this regime of weak coupling, the color-superconduc ting gap (j) can be c omp uted within the fundamental 
' theory of strong interactions, quantum chromodynamics (QCD) IllL ll^ IT^ UJ. It was first noticed in Refs. 
\^ , [Tol IllL IT^i that in order to describe color superconductivity correctly it is crucial to take into account the specific 
• energy and momentum dependence of the gluon propagator in dense quark matter. It turned out that the long-ranged, 
' magnetic gluons generate a logarithmic enhancement in addition to the standard BCS logarithm and by that increase 
-j-jJ ■ the value of the gap at leading logarithmic order. Furthermore, due to the energy and momentum dependence of 
the gluon propagator the gap also is a function of energy and momentum and therefore a complex quantity [l^ ■ 
Complex gap functions are well-known from the investigation of strong-coupling superconductors in condensed matter 
physics for more than 40 years [IE IH, 113, EH . 

The authors of Ref. |l2j | estimated the magnitude of the imaginary part of the color-superconducting gap for 
massless quarks by considering the cut of the magnetic gluon propagator in the complex energy plane. They found 
that Im0 = on the Fermi surface and Imt/) ^ gJiecf) exponentially close to the Fermi surface, |fc — /i| ^ iiexp{~c/ g) 
; ^ • where 5 <C 1 is the QCD coupling constant in the limit of weak coupling. One only has Im ~ Re for quarks farther 
' away from the Fermi surface, |A: — /i| ~ gfi. Therefore, considering quarks exponentially close to the Fermi surface, 
the approximation (j) ~ Ret/) is valid for (f> up to corrections of order g to the prefactor of 0, i.e., up to its subleading 
order. However, the contribution of Im0 to (f> through Ret/) has not been estimated yet. 

In this work we show that, in the 2SC phase and to subleading order, Im0 does not contribute to (j) for quarks 
exponentially close to the Fermi surface. To do so, all momentum and energy dependences of (j) must be included 
in the gap equation. The appropriate starting point for that is an energy- and momentum-dependent Ansatz for 0, 
(/) = (/)(fco, k), where kg and k are treated as independent variables. As known from the solution for Re(/), the integrals 
over energy and momentum in the gap equation yield large logarithms, ln(/i/(/)) ~ which cancel powers of g from 
the quark-gluon vertices. Due to these logarithms one must actually compute or at least carefully estimate these 
integrals in order to determine the importance of the various terms contributing to Hecj) and Im^. Moreover, for a 
complete account, not only the magnetic cut but also the electric cut, as well as the poles of the gluon propagator, 
have to be considered in the solution for lm(f>. In order to illustrate the latter point, for energies just above the gluon 
mass nig ~ g^ one has Im (f> g^^ ^ Re (f> due to a large contribution from the emission of on-shell electric gluons. 
In order to estimate how this term feeds back into Re (f> requires a careful analysis. Treating energy and momentum 
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as independent variables and solving the coupled gap equations for Im0 and Re0 self-consistently is therefore a 
non-trivial problem and, moreover, leads to interesting insights. 

To date, the gap function has never been calculated by treating energy and momentum independently. In Refs. 
[T^ Hsl it was assumed that the off-shell gap function is of the same order of magnitude as the gap on the 
quasiquark mass-shell, 4){kQ,]i) « (/)(efc,k) = (/)k. Furthermore, all contributions that are generated by the energy 
dependence of the gap, i.e., by its non-analyticities along the axis of real energies, are neglected completely against 
the non-analyticities of the gluon propagator and the quark poles. In Ref. |2J| it was pointed out that at least in 
order to calculate corrections of order g to the prefactor of the gap, its off-shell behaviour must be included. In 
Refs. 0,0,0, on the other hand, the gap has been calculated within the Eliashberg theory In this model it 
is assumed that Cooper pairing happens only on the Fermi surface. Following this assumption, the external quark 
momentum in the gap equation is approximated by fc « /i. By additionally assuming isotropy in momentum space, 
the gap becomes completely independent of momentum and a function of energy only, 0(a;, k) « (f)(uj,fi) = 4>{uj). 
Originally, the Eliashberg theory was formulated in order to include retardation effects associated with the phonon 
interaction between electrons in a metal. Since the energy that can be transfered between two electrons by a phonon 
is restricted by the Debye frequency lud, Cooper pairing is restricted to happen only at the Fermi surface |21[ . In 
quark matter, however, such an assumption is problematic if one is interested in understanding quark matter at more 
realistic densities where the coupling between the quarks becomes stronger. In this case, quarks further from the 
Fermi surface also participate in the pairing p5| . If color superconductivity is present in the cores of neutron stars, 
the coupling is certainly strong, (7 ~ 1, and a non-trivial momentum dependence of the gap cannot be neglected. 
The present analysis is, strictly speaking, valid only at weak coupling. However, treating energy and momentum as 
independent variables might still be helpful to catch some aspects of color superconductivity at stronger coupling. 

Besides affecting the value of the energy gap (jj at some order, Im (j) also contributes to the damping of the quasiquark 
excitations in a color superconductor. With respect to the anomalous propagation of quasiquarks, it is shown that 
the imaginary part of the gap broadens the support around the quasiquark poles, see Eq. (|18() below. This broadening 
strengthens the damping due to the imag^inary part of the regular quark self-energy, E, which is present also in the 
non-colorsuperconducting medium "26*, 27, 28, 29, "s^. 

This paper is organized as follows: In Scc^lthc 2SC gap equation is set up within the effective theory derived in Ref. 
[2^. In Sec, mil the gap equation is first decomposed into its real and imaginary parts and then solved to subleading 
order, cf. the schematic outline given after Eq. H33|) . It is shown that Im0 contributes to (j> at sub-subleading order 
for quarks exponentially close to the Fermi surface and that Im0 ~ Ke(f) at |fc — /x| ~ gfj,, which justifies previous 
calculations. Furthermore, an analytical expression for Imcj) is given, see Eq. (|91|l below. In Sec. llVl the conclusions 
and an outlook are given. A somewhat more detailed presentation can be found in Ref. Tel. 

The units are h ~ c — ^ I- 4- vectors are denoted by capital letters, ~ (fco,k), with k being a 3- 
vector of modulus |k| = k and direction k = k/fc. For the summation over Lorentz indices, we employ a metric 
gi^'^ = diag(-|-, — , — , — ) and perform the calculations within a compact Euclidean space-time with volume V/T, where 
V is the 3-volume and T the temperature of the system. Since space-time is compact, energy-momentum space 
is discretized, with sums {T/V)J2k = ^X]n(l/^) X^k- ^ large 3-volume V, the sum over 3-momenta can be 
approximated by an integral, {l/V) J2w — J d'^'^/ {"^tt)^ . For bosons, the sum over n runs over the bosonic Matsubara 
frequencies wjj = 2n7rT, while for fermions, it runs over the fermionic Matsubara frequencies = (2n + 1)ttT. 



II. SETTING UP THE COMPLEX GAP EQUATION 

The complex gap equation is set up within the effective theory derived in Ref. [24j . This has three major advantages 
over a treatment in full QCD: Firstly, self-consistency of the solutions of the Dyson-Schwinger equations for the quark 
and gluon propagators is only required for those momentum modes considered as relevant for the physics of interest. 
In full QCD, on the other hand, self-consistency has to be maintained for all degrees of freedom. Secondly, by a 
special choice of the cutoffs for relevant quarks and gluons, Aq and Agi, one can implement the kinematics of quarks 
scattering along the Fermi surface into the effective theory. Considering quarks with momenta \ k — l i\ < Aq as relevant 
degrees of freedom, one can define the projector onto these modes in Nambu-Gor'kov space as [2J| 

ri{K, Q) ^ ) e(Aq -\k- fc^i) , (1) 

where A|. = (1-1- 6707 •k)/2 projects onto states with positive (e = or negative (e = — ) energy (quark masses being 
neglected). The quark modes far away from the Fermi surface as well as antiquarks have the projector V2 = ^ ~ Vi- 
They are integrated out and are contained in the couplings of the effective theory. For the gluons we introduce the 
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projector 



Qi(Pi,P2) ^ e(Agi-pi)^!.1^, ■ (2) 



Consequently, relevant gluons are those with 3- momenta less than Agi, while gluons with larger momenta, correspond- 
ing to Q2 = I — Qi, are integrated out. Choosing the cutoffs according to 

Aq < gfi < Agi < /i (3) 

the energy of a gluon exchanged between two quarks is restricted by po < Aq. Its momentum, on the other hand, 
can be much larger, since p < Agi. This reflects the fact that quarks typically scatter along the Fermi surface and, 
due to the Pauli principle, do not penetrate deeply into the Fermi sea. In addition to that, gluons with po ^ p have 
the property that they are not screened in the magnetic sector and therefore dominate the interaction among quarks. 
The third advantage of this effective theory is that by expanding the numerous terms in the gap equation in terms of 
Aq/Agi ~ g one can systematically identify contributions of leading, subleading, and sub-sublcading order. This was 
demonstrated explicitely in 24] for the real part of the gap equation. Similarly, also the terms in the complex gap 
equation can be organized in this way. Obviously, the separation of the scales (f) , gfi, and fj, is rigorously valid only 
at asymptotically large values of the quark chemical potential, where g <C 1. In the physically relevant region, /i ^ 1 
GeV and g ~ 1 , th is scale hierarchy breaks down. For that case, more suitable choices for cutoff parameters have 
been suggested |3ll |. 

The Dyson-Schwinger equation for relevant quarks and gluons can be derived in a systematic way using the Cornwall- 
Jackiw-Tomboulis (CJT) formalism 32]. For the quarks one finds 



[G+]-i \ / S+ 
[G-]-i / I $+ 



(4) 



Here [G"'"]""'^ is the inverse tree-level propagator for quarks and [G~]~^ is the corresponding one for charge-conjugate 
quarks. These effective propagators differ from the QCD tree-level propagator [Gq]~^{K) = If ± /ijo by additional 
loops of irrelevant quark and gluon propagators. In Ref. it is shown that to subleading order in the gap equation 
these loops can be neglected, [G"*"]"^ ~ [Gg ]~^. The regular self-energy for (charge-conjugate) quarks is denoted 
as E^. The off-diagonal self-energies <I>^, the gap matrices, connect regular with charge-conjugate quark degrees of 
freedom. A non-zero corresponds to the condensation of quark Cooper pairs. Equation Q can be formally solved 
for^. 



s+ g- 



(5) 



where 



5± = {[G±]-i+S±-$T([GT]-i + ET) '$±1 ^ (6) 
is the propagator describing normal propagation of quasiparticles and their charge- conjugate counterpart, while 

S± =-([GT]-i+ST)"'$±g± (7) 

describes anomalous propagation of quasiparticles, which is possible if the ground state is a color-superconducting 
quark-quark condensate, for details, see Ref. 5|. To subleading order, it is sufficient to approximate the propagator 
of the soft gluons by the HDL-resummed propagator Ahdl instead of solving the corresponding Dyson-Schwinger 
equation 33], while for the hard gluons one may use the free propagator Ao,22 [l^l- The index 22 indicates that this 
propagator describes the propagation of a hard gluon mode. One has in total 

A^a^^) ^ [Ahdl]^: (P) ^(Agi ~p) + [Ao,22]:: (P) eip - Agi) . (8) 
In the mean-field approximation |22j ] the Dyson-Schwinger equation for the gap matrix <I>+(_R') reads 

^g'^J2 ^atiK - Q) i.iT'^r s+(Q) 7.r^ , (9) 

Q 

cf. Eq. (97) in Ref. [23] . As discussed above, in the effective theory the sum runs only over relevant quark momenta, 
/i — Aq<(7</^-|-Aq. Due to the dependence of the gluon propagators Ahdl and A0.22 on the external quark energy 
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momentum K in Eq. (j^jl, the solution 4>(JC)+ must be energy-dependent itself. Hence, solving the gap equation self- 
consistently requires an energy-dependent Ansatz for the gap function. To sublcading order in the gap equation, the 
contribution from the regular self-energies can be subsumed by replacing qo/Z{ko) in the quark propagators 
Isll. where 



-2. 



is the quark wave- function renormalization factor I8(tI | , with 
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Z(fco)= (1 + r In— ) (10) 



3V27r ' 



(11) 



and 



The effect of Im S on Re has been studied in Ref . (29| , where it is shown that Im S suppresses the formation of quark 
Cooper pairs. The corresponding corrections, however, are shown to enter Re cj) only beyond subleading order. In the 
following it will be assumed that Im E enters Re </> through Im (f> also only beyond subleading order. Consequently, 
Im E will be neglected completely. This is self-consistent since it turns out that Im itself contributes only beyond 
subleading order to Hecj). The main contributions to Imt/) are expected to arise from the energy dependence of the 
gluon propagator, and not from Im E. This amounts in neglecting the cut of the logarithm in Eq. H10|l when performing 
the Matsubara sum in the complex gap equation Q. 

For the sake of definiteness, a two-flavor color superconductor is considered, where the color-flavor-spin structure 
of the gap matrix is 

^+{K) = J3T275 A+ e(A, - |fc - (t>{K) . (13) 

The matrices {J?,)ij = —isijs and (T2)/g = —i^fg represent the fact that quark pairs condense in the color-antitriplet, 
flavor-singlet channel. Then the anomalous propagator reads 

S+(Q)- J3r275 Aqe(A,-|q-A^|)^— (14) 
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where 



ek = V(fc - + 02 . (15) 

Here we employed the analytical continuation j^p (jy^ [itL [T9ll2l| . Besides its poles at qo = ±Z{eq)eq = ±iq the 
anomalous propagator S+ obtains further non-analyticities along the real qg-axis through the complex gap function 
(j). As presented more explicitly in Appendix the energy dependence of the gap function (l){K) gives rise to a 
non-trivial spectral density 

P0(w,k) = ^ [0(a'-|-ze,k)-0(w-ie,k)] , (16) 

which is directly related to the imaginary part of the gap function via 

Im 0(0; -)- ie, k) = 7rp0((jj,k). (17) 
For the spectral density of the anomalous quark propagator this yields 

PH(w,q) = ^ [S(tj + iry, q) - S(cj - i?7, q)] 

= ~Z^{u^)V , - sign(c.) Z^(gq) Re 0(c. + »7y, q) 6{u:^ - g^) , (18) 

where the cut of Z{Ld) has been neglected. Also the non-analyticities of can be neglected: In the region \k — /i| ~ 
one has ~ Ret/f, whereas for — /i| ~ ^ it is ~ \q — Hence, in Eq. H18|) and in the gap equation © 
one may write ~ ^^(q — /i)^ + (Re 0)^, which is continuous across the real energy axis. From Eq. (|18|l it becomes 
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obvious that p^{u},ci) ^ leads to a broadening of pH(u;,q) around the quasiparticle pole at w = e^. However, in 
order to describe the damping of quasiquarks self-consistently, it is necessary to include ImE. Further interesting 
details of the damping due to Im can be found in the recent analysis of Ref. 

Inserting Eq. H14|l into Eq. multiplying from both sides with J3T275Aj^, and tracing over color, flavor, and Dirac 
degrees of freedom, one finds with [Aq, 22106 = 



Sab A^^22 and [AHDLlafc = ^^b ^'huL 



'^(^) = T y E Ti-^ {Ki,K^^) ^''"{K - Q) HQ) HQ) , 

Q 



where the remaining traces run over the internal Dirac indices. For convenience, the quark propagator 

9o - [Zilo) <^qY 2eq ^ qo - crZ(go) Cq 



A(g) 



(19) 



(20) 



has been introduced. Before the actual solution of Eq. (|19|l in the next section, we first discuss the power-counting 
scheme of the gap equation in weak coupling. In order to fulfill the equality in Eq. (|19|) , the integration over energy 
and momentum on the r.h.s. must yield terms of the order (j^/g^. After combination with the prefactor they are of 
order cf), which is the leading order in the gap equation. Accordingly, terms of order g cf) are of subleading order and 
terms of order g'^cj) are of sub-subleading order. Until now, the color-superconducting gap is known up to subleading 
order, i.e., up to corrections of order g to the prefactor of the gap 0. 

III. SOLVING THE COMPLEX GAP EQUATION 

A. Derivation of the coupled gap equations of Re (j> and Im (j> 

In order to derive of the coupled gap equations of Re (p and Im <j) we first rewrite the Matsubara sum in Eq. (|19|l as 
a contour integral, 

(21) 



X^^*(fco,P,q)^r A'^\Q-K)K{Q)m)= I ^^tanh(^) A^^*(Q-if)A(g)0(Q) 



where p = q — k and A^ denotes the longitudinal and A' the transverse gluon propagator in pure Coulomb gauge, 
cf. Add. IbI The contribution at qa — ko, where the cut of the gluon propagator is located, has to be omitted. The 
contour C is shown in cf. Fig. ^ In order to introduce the spectral densities of the gap function and of the magnetic 
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FIG. 1: The contour C in Eq. 12111 encloses the poles of tanh[go/(2r)] on the imaginary go axis. The additional poles and the 
cut at qo = ko arise from the gluon propagator, while the two poles on the real axis are due to the quasiquarks, cf. Eq. 1181 . 
The yet undetermined non-analyticities of the gap function on the real go— axis are indicated by the shaded area. 



and longitudinal gluon propagators the contour C is deformed corresponding to Fig. El The contour integral can now 



FIG. 2: Deforming the contour C to introduce the spectral densities of A^"' and cj>, cf. Eq. 12211 . 



be decomposed into three parts 



M^^*(fco,P,q)=/4* + /o"+4*. (22) 



The first part, loo, is the integral along a circle of asymptotically large radius. It can be estimated after parameterizing 
dqo — i\qo\e^^dd and considering the limit |go| ~> oo- To this end we write without loss of generality, cf. Appendix 1X1 

^{K) = Uk) + 0(k) , (23) 

so that the energy dependence of (t>{K) is contained in (l){K). For asymptotically large energies, |A:o| oo, we require 
that (i){K) and (p{K) <^(k), where (j){'k) is a real function of k only. Furthermore, we have A^{Q~K) |q— k|~^ 
and A*(Q — K) —Qo'^, cf. Eq. IIB2p . It follows that the longitudinal contribution ~ q^^ and the transversal 
^oo ~ ^-iid hence that in total — > 0. The integral /q'* runs around the axis of real qq, 

oc 

l"= J ^itanh(^) A^^*(Q-if) [A(go + i?/,q)</'(<7o + *??,q)-A(<Zo-iry,q)0(go-i?7,q)] ■ (24) 
Applying the Dirac identity and using Eq. (jlGf) and (j){uj + if]) + (j){uj — irf) — 2Ke(j)(ijj) one obtains 

oo 

= ^E-^-. / rf'/oJtanh(|i)A"(,o-fco,p)^^(<Zo)^^^ 
2eq J 2 \2TJ qo-creq 



^^tanh (^^j Z2(?q)Re0(e~q,q) ^ A'-\ai^ - fco,p) , (25) 



where Vx denotes the principal value with respect to the pole at x. The first term arises from the non-analyticities of 
the gap function, cf. Eq. and the second from the poles of the quark propagator at Qq = ±eq, cf. Eq. (|2nj|. 

The last integral I^'^ circumvents the non-analyticities of the gluon propagator as well as the pole of tanh. One 
finds 

oc 

-^J^aT'o / d<Zoicoth(|^)</.(qo + fco,q)^^(go + M^-^%^ 
2eq ^ J 2 \2TJ qo-creq + ko 

— OO 

-TA^'*(O,p)A(fco,q)0(fco,q) . (26) 

The first term is due to the non-analyticities of the longitudinal and magnetic gluon propagators, A^'*. The second 
arises from the pole of coth(go/2r) at go = corresponds to the large occupation number density of gluons in 
the classical limit, q^ <C T. This contribution has been shown to be beyond subleading order T^l and will therefore 
be discarded in the following. 



'fco 

^ CT=± 
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After analytical continuation, fco — > uj + ir], the Dirac identity is employed in order to split the complex gap equation 
(|19|1 into its real and imaginary part. For its imaginary part one finds using Eqs. 125I2()|) and using the fact that 
Re0((7o) and Z^{qo) are even functions 



Im7W^'*(tJ + M7,p,q) = -^Re0(?q,q)Z2(?q) ^ cr/^*(t^-CTeq,p) 



tanh 



2T 



coth 



— > Cr-P 



1 r,= 



p^'\uj ~ go,p)p0(go,q) 



ago : ^ (go) 

f-t/ , ■ \ , T KAl.tl 



tanh (|;) 



coth 



2T 

^ - go 

2T 



ImA^_^ + IT], p, q) + ImAlg (w + i?/, p, q) 



(27) 



The first term on the r.h.s of Eq. H27|) . ImAI^*, contains all contributions to (\>{K') that have already been considered 
in previous solutions to subleading order. Here the gap function is always on the quasiparticle mass-shell. The second 
term, ImAlg*, is due to the non-analyticities of (f){K). It has been neglected in all previous treatments due to the 
approximation 



/3H(t^, q) ^ -sign(cj) Re 0(eq + i?7, q) Z'^{eq)5{, 



(28) 



cf. Eq. (41) in Ref. By doing so, one always forces the gap function on the r.h.s. of the gap equation on the 

quasiparticle mass-shell. The gap equation takes then the standard form, cf. Eq. (19) of Ref. 34] . The occurrence of 
the external energy ik on the r.h.s. due to the energy-dependent gluon propagators indicates that the solution still 
possesses some energy dependence although not provided by the Ansatz. 

In the limit of small temperatures, T — > 0, the hyperbolic functions in Eq. H27|l simplify, yielding for w > 



Im7W^*^(,(tj + ir/, p, q) 



TT 



Z2(e~q) Re0(?q, q) p'^^LO - ?q, p) - ?q) 

/'^(cj - go,p)p0(go,q) 



^'(go 



= ImA^^*y^g(a; + i??, p, q) + ImA^^*j,^Q(Lj + i-q, p, q) . 



Inserting the Matsubara sums IvaMj^ g back into Eq. ifT^I) . 



a f da 

Im(/)(tj + M7,k) = YJ J^Y^^^'^^'^'P'I^ [Im7W^(w + M7,p,q) +ImA1e(cj + M/,p,q)] 

= AiiJ + iri, k) -f B{uj + i7], k) , 
we have finally derived the gap equation for lni(j){Lu + irj, k). The traces over Dirac space are given by 



Tr5(fc,P,g) = 



Tr'(fc,_p, g) 



{k + qf 



2kq 
2k q 2k qp^ 



(29) 
(30) 

(31) 
(32) 

(33a) 
(33b) 



We found that the imaginary part of the gap function can be split into a term A{u! + irj, k) which contains all 
contributions of the real part of the gap function on the quasiquark mass-shell, and a term B{uj + ir], k) which contains 
all new contributions with the imaginary part of the gap function off the quasiparticle mass-shell, cf . Eq. (|32|l . In the 
following it will be checked whether the known solution for Re0(ek, k), which neglects the contributions contained in 
B, is self-consistent to subleading order. To this end, we can use the leading order solution for Re(/)(ek,k), which is 
given by 



<j){y) = (j) sin l^— j , 

where (j) denotes the value of the gap function on the Fermi surface to leading logarithmic order in o [Tc 



Mexpl-- 



(34) 



(35) 
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^ Ai 


A 

UJ ~ Al>j;>g 


, . A 

U ~ i\-§>j;>0 


uj ~ rrig + Al>y>0 


^; r:; 


dominant gluons 


t-cut 


t-cut 


t, ^-cut 


^-pole 


t-pole 


A 




3 COS (^) 






g<t> 


B 


g'A 


g'A 




g^A 


gA 


H[A\ 








<i> 


g<t> 


n[B] 




g^c^ 




g^4> 





TABLE I: Estimates for the terms A and B, cf. Eq. lEU, and for Tl! [A] and [S], cf. Eq. (E7J, for different energy scales and 
|fc — ^1 -C M . The gluon sectors dominating the respective energies are indicated. 



The variable < y < 1 defines the distance from the Fermi surface through the mixed scale 

Ky = (f^yM^-y , (36) 

where M ~ g/x is defined in Eq. (|12|l . A given value of y corresponds to momenta k with |fc — /i| ^ Aj,. Note 
that Ki = (j) and Aq = M. Furthermore, we have Kg ^ e~'^^^ M, which is smaller but still of the order of M. 
Correspondingly, we refer to quarks with momenta \k — fi\ ^ Ai<y<g as exponentially close to the Fermi surface and 
to quarks with |A: — /x| ~ Ag>y>o as farther away from the Fermi surface. Inserting (j){y) into A{uj + i??, k) one can 
estimate A{uj + ir], k) for different energy and momentum regimes. This is done in Sec. IIII Bl In the second iteration 
the part B{oj + irj, k) is estimated by inserting ~ A/tt into the expression for B{lu + irj, k), which is done in Sec. 
IIII CI In Sec. IIII D1 these estimates are used to write 



Re0(ek,k) = -V 

IT 



As Ao 



A. 



y -^(^+^'^^k) + g(c.+.,,k) 

(T = ± 



cf. Eq. l|A7a|l . where the integral over uj has been split according to the different energy regimes of the estimates 
for A and B, cf. Table Then, according to the discussion after Eq. (|20ll . terms of order g"-(f> contribute to the 
(sub)"-leading order to Re0(ek,k). The main results of this analysis are summarized in Table^for momenta close 
to the Fermi surface, |fc — <C M . The columns correspond to the various energy regimes of these estimates. In the 
first line the dominant gluon sectors are given. The cut of the transversal gluons gives the dominant contribution to 
A and B for energies smaller than the scale M. At the scale M, the longitudinal and the transversal cut contribute 
with the same magnitude. At that energy scale also the poles of the gluons start to contribute as soon as w > nig. 
At energies just above rUg the longitudinal pole dominates over the magnetic pole, whereas for larger energies up to 
2/i the transversal pole gives the leading contribution. The respective orders of magnitude of A and B, estimated at 
the various energy scales, are given in the subsequent rows. It is found that either B ~ gA or i3 ~ g^A and therefore 
B <s:^ A. Finally, the orders of magnitudes of the contributions to the Hilbert transforms Ti [A] and Ti [B] , obtained by 
integrating over the respective energy scales, cf. Eq. (|37|l . are listed. When calculating the Hilbert transform H [A] 
up to energies uj ~ Ag, the transversal cut gives a contribution of order 0, which is of leading order. Since B ^ g^A 
for these energies, Ti [B] ~ g'^cj), i.e., the contributions from B are beyond subleading order. Integrating over larger 
energies uj ~ M, it turns out that Ti. [A] gcf), which gives a subleading-order contribution. Since B ~ gA in this 
energy regime, the corresponding contribution from Ti [B\ is again beyond subleading order. For the contributions 
from the poles one finds that, in the region oj = vig + Ai>y>o, a contribution of order cf) is generated by Ti [A\ (which 
combines with (?!)(k) to give a contribution of order g4> \n total, i.e., of subleading order). Since in this energy regime 
B ^ g^A, the corresponding contribution Ti, [B] is again only of order g^(j>, i.e., beyond subleading order. Integrating 
over large energies up to 2/i, Ti. [A] gives a contribution of order g (jj, which is of subleading order. Since in this regime 
B ^ gA, it follows that H [B] is beyond subleading order. 

In Sec. IIII El it is shown that also (/)(k) ~ cj) and that the imaginary part of the gap function again contributes 
to 0(k) only at 5^0, i.e., beyond subleading order. This finally proves that the imaginary part of the gap function 
enters the real part only beyond subleading order. Hence, to subleading accuracy, the real part of the gap equation 
can be solved self-consistently neglecting the imaginary part of the gap function for quark momenta exponentially 
close to the Fermi surface. On the other hand, the real part of the gap function enters the imaginary part of the 
gap function always to leading order. For energies, for which B ~ g^A, the imaginary part of the gap function can 
be calculated to subleading order from the real part alone. In particular, this is the case for the regime of small 
energies. Furthermore, since for these energies only the cut of magnetic gluons contributes, Imc/) can be calculated to 
subleading order without much effort, which is done in Sec. IIII d In Sec. IIII Fl we reproduce the known gap equation 
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for Re(/()(ek,k) by Hilbert transforming A{lli + 177, k), cf. Eq. (jA7a|) . and adding the energy independent gap function 
4){k), of. Eq. 1231. 

B. Estimating A 

For the purpose of estimating the various terms contributing to A, cf. Eqs. (|29I32(I . one may restrict oneself to the 
leading contribution of the Dirac traces in Eq. H19|l . which is of order one. The integral over the absolute magnitude 
of the quark momentum is / dqq^, while the angular integration is / dcos9 = J dpp/{kq). Furthermore, we estimate 
Z^(eq) ~ 1. The contribution to A from M^_^t=o 123) > which arises from the cut of the soft gluon propagator, 

is 

S Agi 

XutC'^.k) ~ .g^y ^Re0(eq,q) J dp p p^^^^{uj* , p) , (38) 

A 

where lo* = lo — and < lu* < uj. Furthermore, S = min(ti;, Aq) and A = max(|^ — with = k — fi and 

^ = q — fjL. In the effective theory, both |^| and |^| are bounded by Aq ~ gjjL. Due to the condition \ < p < Agi it 
follows that ^cut = for u) > Agi + Aq ~ /i. For the purpose of power counting it is sufficient to employ the following 
approximative forms for cf. Eq. HB5blB5d|l . 

,p) (39a) 

pUc.*,p) J . (39b) 

These approximations reproduce the correct behavior for w* <C p ^ Wg, while for uj* < p < Agi they give at least the 
right order of magnitude. With that the integration over p can be performed analytically. For energies ui < Agi one 
finds for the transverse part 



J Cq 





arctan — ^ — — arctan 



(40) 



For all C ^ and to < Agi one has Agj/Af^w* ^ 1 and the first arctangent in the squared brackets may be set equal 
to 7r/2. In the case that lu ^ M one has \^ /NPuj* 3> 1 and the two arctangents cancel. If w ~ Af, on the other hand, 
the arctangents combine to a number of order 1. Finally, in the case that lo ^ Ky with < 2/ < 1, one has always 
< ^ < Ay due to the theta- function in Eq. H29|) . Moreover, it turns out that the arctangents cancel if C > ^y/Si 
since then X'^/iM'^io*) ~ C/[M^{uj - C)] > 1 for all ^<uj. 

In the case that oj ^ (p, the integral over ^ does not yield the BCS logarithm, cf. Appendix lO and we find 

^*„t(0,k) ~ 52^. (41) 

For larger energies uj ^ Ay with < y < 1 one substitutes £,{y') = Ayr, d^/^ — ln(0/A'/) dy' and obtains with Eq. H34f) 



V 

^cut(^^>k) - \n (^j-^ (j) j dy sva{^^ ^ gcj) cos 



(42) 



Here, the integration over < ^ < Ai has been neglected, since it gives at most a contribution of order g^^, cf. Eq. 
In the longitudinal sector one finds for the integral over the gluon momentum p 

1 A2 - AAgi + X'^ 




1 










arctan 




— arctan 











Ah A2 + X2 



, for A < X 
, for A > X 



(43) 
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UJ ^ (j) 




a; ~ M 


M <^uj < ^ 






3 COS (2^) 


g<t> 











g(t> 





TABLE 11: Estimates for ^^ut difTerent energy scales and C, <^ M. 





LU ^ (f) 




u ^ M 


M <^uj < ^ 


•^cut 








g<P 





A^ . 






g<t> 


9<l> 





TABLE IIL Estimates for A''^^^ at different energy scales and C, ^ M. 



where = Sm^. Since C, < Aq ^ X, solely the magnitude of cu decides whether A < X or A ^ X is realized. 
It follows that, in contrast to the transversal case, the order of magnitude of A^^^ is independent of Energies 
Lu ~ Ao<j,<i correspond to A < X. For the special case w ~ Ai one finds 



^cut(w,k) 



-Re</>(e„q)-^g 0-. 



(44) 



For Lo ^ Ao<y<i one obtains 



-4cut(w,k) 



dy' sin ( ^ ) ^ 9 (1^ I dy' sin ( ^ 



Try 



M 



M 
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(45) 



Hence, for y > the term A^^^^ is suppressed by a factor {(t)/M)y as compared to -4*ut whereas for j/ = the 
longitudinal and the transversal cut contribute at the same order, -4^^^ ^ -^cut ^ 9 'P- For much larger energies, 
M <C w < fi, we have A = w* ~ w ^ X (note that C is bounded by Aq). It follows with J = Aq ^ Ag 



u 

Ai^d^,k) ^ g<pj dy sin ' J^i 



1-- 



(46) 



Hence, A^^t ^ ■^lut this large-energy regime. The results for A^'^^^ are summarized in Tables IhI and IIIII 
The contributions from the gluon poles to A read analogously to Eq. (|38|l 

S 2)1 

4'l(a;,k) ~ g2y'^Re0(eq,q) j dpp/pl(c.*, p) - c.„(p)] . (47) 

l?-CI 

The boundary p < 2fj. in the integral over p is due to the constraint ^ < Aq, cf. Fig. 13 Therefore, also w^.f < 2/i. It 
follows that .4p*[^ = for uj > since for these energies one has w* ~ w > LUg^t always and the 5— function in Eq. 
(|47|l vanishes. Hence, we can restrict the analysis of Eq. (|47ll to nig < uj < 2/i. 
For the transverse sector we approximate 



Ppolo('^*(P)'P) - - 



and uJt (p) — \ p^ + for all momenta p and obtain 



2wt(p) 



2p 



(48) 



4„i,(l^, k) - 9^ jf- Re 0(eq, q) j dtot S{u;* - uot) . 



(49) 
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FIG. 3: Hard gluon exchange with momentum p > Agi ~ fi. The quark has to remain within a layer of width 2Aq ~ p/u around 
the Fermi surface. This effectively restricts the hard gluon momentum to p < 2/i + 2Aq < 2fj,. 



The condition uj* = LUf can be satisfied only if w > y + C^. Furthermore, the condition uj* = ujt sets an upper 
boundary to the integral over ^ given by ^max = niin{Aq, (w^ — — C^)/2(a; — ()}. Hence, the BCS logarithm is 



generated for energies u ~ + + with < y < 1, since then ^max ~ A^. For such energies we have 



Jma,x 

4oie('^.k) ~ i?' y ^ Re 0(eq,q) cos (^^ 



(50) 



For larger energies up to 2/i one has Apgi^{cj, \s.) g^. 
In the longitudinal gluon sector we approximate 

pU('-Kp),p)--^ (51) 

for gluon momenta p not much larger than nig. Such values of p are guaranteed if we consider energies of the form 
LO — ^Jrrig + + Ayj with ^ < yi < 1. As in the transversal case we simplify ^^(p) ~ ^p"^ + nig and find 

4„i,(a;,k) ~ r ^Re^(^^^q) / dco, , ^ 6{u* - toi) , (52) 

where ^max is defined as in the transversal case. For the considered energies one has ^max ~ A^^ . We find 

2 f '^^^^Ai. „\ ('^-^q) 



-^poie(^.k) ~ 5^ / — Re^!>(eq,q; 







(w - eq)2 - mi 



w — nio 



cos 



(^) . (^^) 



where ^max and w > to^ was exploited in order to simplify the fraction under the integral, a;^/(a;^ — m^) ~ 

uu/{uj — nig). Furthermore, the BCS logarithm has cancelled one power of g. Hence, for Q ~ nig one has A^^^^^ ^ g(f>- 
For quarks exponentially close to the Fermi surface with ( ^ ^y2/2 and 2/2 > 2 5, we find 

'My 



Al,M k) ~ 5 (^) ' COS (^) , (54) 
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^ < \Jm1 + 






'^polc 







g<t> 


A^ , 





5 0(^)''cos(^) 





TABLE IV: Estimates for A^^^^ at different energy scales and for ~ A^^, where < j/2 < 1 and y = min{j/i , 2/2}- 



where y = minlyi, 2/2}- For much larger energies, lo 3> Trig, we have to consider gluon momenta p 3> mg, for which 
the spectral density of longitudinal gluons is exponentially suppressed 



Ppolc(^^i(p),P) 



p 



We find for nig <C ijJ < 2/i with uji{p) ~ p 



(55) 



Aq ^ 



(5(ti;* - Lj^) 



Ao 

9 (P J Y 



2(c.-02n 



3m2 



ff 9> exp 



2^2 
3m2 



(56) 



which is the continuation of Eq. (|53|l to large energies and for all C ^ ^q- The results for A^^^^ are summarized in 
Table HVl In the following Sec. IIIICJl these results will be inserted into Eq. for in order to estimate B. 



C. Estimating B 



(7A{qQ,q) 



90 



(57) 



The term B in Eq. H32|) can be estimated using 

Scut(w,k) ^ — dqoY^ 

'^-^ A 

for the Landau-damped gluon sector. We introduced to' = uj ^ qo < uj and A = max(|^ — (^\,lu'). Furthermore, we set 
Z'^{u!) ~ 1. From the condition A < Agi in Eq. lf57|) and from A{qo, q) = for qo > 2/i it follows that K^^\(w, k) = 
for lli > Agi + 2fj, ^ 3p,. Inserting the approximative forms (|39|l for into Eq. (|57|l the integration over p can be 
performed analogously to Eqs. H40I43() . In the transverse case one finds 



^ [dqoY. 



(TA{qo,q) 



arctan 



APu. 



arctan 



A3 



(58) 



Analogously to Eq. H4()|l . one first determines the domains of w, Cj and ^, where the arctangents in the squared brackets 
do not cancel. Since w' < < 3/x it is Agj/M^o;' ^ 1 and the first arctangent in the squared brackets may be set 
equal to 7r/2. Furthermore, one finds that the argument of the second arctangent is not very large as long as the 
conditions w — M ^ go and |^ — CP ^ M'^ioj — qo) are fulfilled. In order to satisfy the first condition, we restrict the 
integral over go to the region max{0, uj — M} < q^ < lo. The second condition becomes less restrictive if simplified to 
1^ — < Af^w. The resulting estimate for will turn out to be small so that a more elaborate estimate is not 
necessary. For energies oj ^ 4> one may use Eq. (|41|l to estimate A ~ ^'ut ~ 9'^4>- For ( <^ M one has 



<>,k)~5' 



In 



5V 



(59) 
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The logarithm under the integral prevents the generation of the BCS logarithm, cf. Appendix |0 For C < M the 
integration over ^ is restricted to the region |^ — CI < ^i/?,- This yields the estimate k) ^ (A1/3/M) ~ 

For Lo ^ Ky with < y < 1 we conservatively estimate A ^ -^cut ^ 9't'i cf . Eq. 142|) , and obtain similarly to Eq. 



e*,t(A„k)^g3^ / ^In 



^ ~ A 



eq + Ay 



(60) 



where we assumed C ^ M. Again no BCS logarithm was generated due to the additional logarithm. For C, M the 
integration over ^ is restricted to the region |C — CI < If conservatively estimate A ^ g cj) throughout this 

region, we obtain k) ~ cj) {(j) / M)y / ^ . 

For energies lu ^ M and larger, the condition |^ — CI'"* ^ M^w is fulfilled for all ^ < Aq. Considering ^ nig + Ay 
with < y < 1 we find that the dominant contribution comes from A ^ -^p^j^, cf. Eq. (|54|l . when integrating over 

Vrig + Kl < qO < mg + Ay 



m„+A„ 



Kutimg + Ay,k)^g^(t)[^ [ dqoY] 

n ™ I A o — in 



9o 



mg+Ai 



90 - rUg 
Ao 



5> 



E 



^2 _ ^2 



(61) 



The contributions from the other gluon sectors are estimated with A ^ g (j)^ cf. Eqs. (|46I52|1 . 



/ ^ / dqo y: — I 

J <^ci J ^ go - creq J 

uj-M "^-^^ 



g 



eq 

MV 



In 



-In 



UJ- M + Ea 



LO — M — e„ 



(62) 



where the logarithms again prevent the generation of the BCS logarithm. The factor (M/w)^ arises from expanding 
the logarithms for uj 3> M. Hence, for energies of the form uj ^ nig + Aj, with < J/ < 1 we have 



while for larger energies A^^^-^^ does not contribute anymore and S*^^(aj, k) ~ g^(j) [MjuSf' , cf. Eq. H62|) . 
For the cut of the longitudinal gluons one obtains 



(63) 



eut(^,k) 
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^q ^-j- 



7~± 9o-f^eq 

where 1(A) is defined in Eq. H43|l . Analogously to the analysis of y^^ut one finds ior uj ^ (f> 



(64) 



,k) 



— dqo y 77 3 T7 / ~ 



/0 



M 



and similarly for ~ Aj, with < y < 1 



(65) 



S^„t(A„k)~g^</.(^ 



(66) 
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Ul ^ (f> 




to ^ rUg + Al>y>0 












9'M^r 













TABLE V: Estimates for B^;; at different energy scales and ( <^ M. 





UJ ^ (f> 


UJ ~ Al>y>0 


LO ^ nig + Al>y>0 






9'<i>m''' 








'-'cut 


3 '?!'X7 




/0 





TABLE VL Estimates for S^;,; at different energy scales and C, M . 



For bj ^ rag + Aj, with Q < y < 1 we can simplify w'X(A) ~ 1 and find as in the transversal case, cf. Eq. H61() . 

SLK+A„k) ^ 52^. (67) 

In the limit of large energies, M lj ^ Agi ^ /i, wc estimate the integral over the range u) — Agi < < uo hy 
substituting A ^ g4> and obtain with uj' 2{\) ^ 1 



J J qo - creq J J 



dqo 



MAa 



M 



(68) 



Hence, also B";,^^ becomes small in the limit of large energies. The estimates for S^^^ are summarized in Tables IVl and 

ED 

In the undamped gluon sector the term AlgV=o P^^l* gives the contribution 



ilq 

6^'oU-,k)^.g^/^/d<ZoE 







gp0(go,q) 



2/j 



(69) 



Due to the restriction -p < 2fi it follows with similar arguments as for B^'^^^ that B^^i^{uj, k) = for cj > 4/i. For the 
transversal sector we employ analogous approximations as for .4* j^, and obtain 



A, 



aA{qQ,q) 



2m 



dWf S{ll>' ~ LOt) 



(70) 



This contribution is non-zero only if w > mg. First we consider energies lo ^ nig + A2y-^ and C. ^ ^V2i and analyze the 
two cases yi < y2 and yi > y2 separately. In the first case, the condition uj > ^ + |C ~ CP requires < ^ < Aj,j^ 
and consequently 



S*,i,(a;,k) ^ g^d, I ^ In 



9 



(71) 



where ^ ~ g(/> and the logarithm prevents the BCS logarithm. The second case, yi > j/2, leads to the condition 
Aj,2 — Ay^ < S, < Aj,2 + Aj,j and we find 



A„,+A„ 



'-^ h. 



mo - i/m^ + A^^ + Ay^ 



VI 



(72) 
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where in the last step the logarithm was estimated to be of order {(fi/M)^^ and the integral over ^ to be of order 

For to > 2TOg the upper boundary for the integral over ^ is given by Aq without further restrictions. The integral 
over qo runs over values qo > mg and therefore receives contributions from A^^^^, cf. Eq. (|54|) . As a consequence one 
finds 

Bl,U^,k)^g^^ (73) 

analogously to Eq. H61() . For uj ^ irig the additional contributions from 2mq < qo < u are only ^ g^<p, as can be seen 
in the same way as in Eq. Ht)2|) . and therefore B^^^^ ^ 0^4). However, for > 2/i + 2mg the condition lu' = uo^ can be 
fulfilled only for q^ > lo — 2pLf!^ 2mg, where g(t>, and one finds 



Aq lO Aq uj 



dqo 



Ai Lj — 2fj, Ai 

g30^1^g30_. (74) 



UJ 



In the longitudinal sector the analysis starts similarly with energies uj ^ nig + A2j,i and ^ ~ Ky^ . Since this restricts 
the gluon momentum to p ^ m^, we apply Eq. H51() and obtain 



In the case that y\ < y2, the condition uj > m?g + |^ — Cp requires < ^ < A^^ and one finds 



Ai 



aA{qo,q) (w - go) 



2 



(76) 



Since !i ^ ^ y "^s ^ 1^ ^ CP < ^2j/i is much smaller than ui ^ mg + A2j,i, we neglect qo against uj on the r.h.s of 
Eq. H76|l . Furthermore, we estimate A ^ g (j) and obtain similarly to Eq. I|7HI 

BU^{uj,k) ^ g^^^^^g^^ (^)''^ . (77) 

u ~ nig \ <P J 

The case that yi > y2 leads to the condition Ay^ — A^^ < £, < Aj,^ + A^^ , and we find similarly to Eq. (|72|l 

For larger energies uj ^ 2mg the upper boundary of the integral over qo will just exceed nig, where it is ^ ~ -^poie' 
cf. Eq. One finds analogously to Sp^jg, cf. Eq. ifT^ . that this gives the main contribution of order 

Bl,,,iLj,k) ^ g'cj, . (79) 

In order to estimate B^^^^ for energies uj > 2fi + 2m g, for which 2mg < qo < uj and A is only ^ g (p, we have to employ 
Eq. ^ 

Aq LU 2/J. 

Bl,^{cj,k) ^ g^f^ [ dqo E / dc., exp f-|4) '^C-' - o;,) 

2m„ "^-^^ ™ 



Ao 



9^jd^ j ^ exp 

Ai 2mg 



2(^- go) 
3m2 



M 



2 



5^0 — . (80) 



In the last step, the integral over qo was restricted to the region uj — mg ^ qo < uj due to the exponential function. 
The estimates for B^^^^ are summarized in Tables IVIll and IVIIII 

In the following the estimates for A and B are used to determine the order of magnitude of T-C[A] and Ti.[B]. 
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LO < TUg + Al 


UJ ^ irig + Al>y>0 


2mg < UJ < 2fi 


2fi < UJ < 4fi 


*-^polc 











*-^polc 







g'cb 


9'^ i^f 



TABLE VII; Estimates for Bp'^'j^ at diflterent energy scales and C, <^ M. 





UJ < mg + k\ 




2mg < u! < 2fj, 


2fi < UJ < 4fi 


^polc 





g <i> (tt) 


g'(^ 




R* 

^polc 







3^0 


9'<t> {^f 



TABLE VIII: Estimates for 5p'^*i^ at different energy scales and C, < M . 



D. Estimating HiA] and T-L[B] 

In order to determine the order of magnitude of Re </> and the order of the corrections due to B we estimate the 
Hilbert transforms T-l[A\ and T-l[B]. For that the quark momentum has to be exponentially close to the Fermi surface, 
C <C M , because for quarks farther away from the Fermi surface Im (j) cannot be treated as a correction anymore. 
Furthermore, in that case the normal self-energy E has to be accounted for self-consistently, which is beyond the scope 
of this work. The integral over uj in the Hilbert transforms TL[A\ and 'H[B] is split into the energy regimes which were 
used to estimate A and B, cf. Eq. IpTzjl . As explained in the discussion after Eq. IpTz)! we select for each energy regime 
the most dominant gluon sectors and estimate their respective contributions to Re and Re 0. 

At the smallest scale, < w < Ai, one has lm(f) ^ -^cut ^ 9'^4>i cf- Eq. (|41|l . which yields 

r1d.Y.^^^^^-9^^A^)-9^^. (81) 
J ^ uj-aei, \ekj 

At the scale Ai < a; < Ag one has Im0 ^ -^cya ^ 9 4>: cf. Eq. (|42|l . With the substitution duj/uj = ln(0/A/) dy one 
finds 

As As g 



Al '^=± Al ^1 

The contribution from yl^ut ~ 9 4'{4'/M)y at the same scale, cf. Eq. H45(l . can be shown to be much smaller, 

Al ''-^ 1 

At the scale Ag < w < Aq one has Im0 ^ -^cut ^ ^cut ^ 9 4'i cf. Eq. H42II45|) . and finds 

V td. ^ ^^^^ ^gJ^^Jdy^g^. (84) 



, U! — creic J IjJ 



For energies uj ^ rug + Ay with < y < 1 one has lm(j> ~ -^p^j^, ^ g <j> {M/cj))^ , cf. Eq. H54|l . and one finds with 
doj = ln(0/Af) Ay dy 



V 



mg+Ai 
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For the regime nig < w < 2/i we have Im (j) ~ ^poio ^ 9 4'^ cf. Eq. H50() , and obtain 



™ cr— ± ^ 



Finally, integrating over 2fi < uj < ifi with Imcj) ~ Spoio 5'^'/' {M/oj), of. Eqs. 168I74|I . one obtains 



From Eqs. (|82() and H85|l we conclude that Re0 ~ 0. Furthermore, 7i[S] contributes to Re0 only at sub-subleading 
order. The corresponding corrections arise from the following sources. The first is .B*^^, which is ^ g'^A*^^^ for u) ^ Ky 
with g < y < 1. After Hilbert transformation it yields a contribution of order g'^cf) to Re0, cf. Eq. (|82|l . and is 
therefore of sub-subleading order. For a; ^ Aj^ with < y < g we have ySf^j ~ 5-4fut- From Eq. (|84|l it follows that 
l-L[Bl^^ and are of sub-subleading order. For w = nig + Ky with < y < one has B^^^y^ ^ ''^^^ ^ 

sub-subleading-order contribution seems possible, since the corresponding contribution from A^^^i^ is ^ (j), cf. 1)85(1 . As 
the latter, however, combines with (j) to a. subleading order term, cf. Sec. 1111 V\ it would be interesting to investigate 
if also B^^Yc finds an analogous partner to cancel similarly. Moreover, we found that Sp^j^, ^ gA^^^^^ for nig < lo < 2/i. 
From the estimate in Eq. H8t)|) we conclude that the corresponding contributions to Re <j) are of sub-subleading order. 
The results are summarized in Table ^ In the next section it is analyzed at which order Im (j) contributes to the local 
part of the gap function, 0. 



E. The contribution of Imr/> to 



The gap equation for the energy-independent part 0(k) is obtained by considering the integrals Jq and 1^^, cf. Eqs. 
(|25I26|I . in the limit |fco| oo. Since p ^ 2/i, the gluon spectral densities p^'*(9oiP) are nonzero only for go ^ 2/i. 
Consequently, the integral over in Eq. (|26|l is bounded by —2// < go < 2/i. Then, due to the energy denominator 
under the integral, I^^ tends to zero as l/|fco| for |fco| oo. In the second term on the r.h.s. of Eq. I|25|) one has 
Eq < Aq ~ g/i. It follows that for ko ^ 2fi > p the transverse gluon propagator becomes A* ~ l/fco- In the 
longitudinal sector one has —l/p^. Hence, in the limit |fco| oo only the longitudinal contribution of the 

considered term does not vanish. Smilarly, one can argue that also in the first term on the r.h.s. of Eq. H25|) only the 
contribution from the static electric gluon propagator remains. Consequently, we find for (/)(k) 

3(27r)^ J eq J P 



Re</)(e~q,q) Z^(e~q) tanh ( ^ 



l?-CI 



Z2(w)tanh 



\2TJ 



(88) 



In the limit T —^ the hyperbolic functions simplify. After performing the integral over p and with Tr^(A;,p, g) ~ 1 
and Z^{u!) ^ 1 we obtain 



2/i 



le-ci 



oo 

Re0(eq,q)+7' J 



cr=± 



■ UJ 



(89) 



where the large logarithm arises from the p— integral. With that and assuming C ^ A/, the integral containing 
Re0(eq,q) is found to be of order (p, and hence 0(k) ^ <j). The remaining contribution from is identical to Eq. 
((37|l up to the extra a due to the hyperbolic tangent. One can conservatively estimate this term by approximating 
~ 5 (/) for < < 4/i and all Ai < C < ^o, and adding ^ g cf) [M/ (pY in the range u ~ nig -I- Ay , 1 > y > 0. 
We find that the contributions from /j^ to 4> are of order g'^C) 
proof that the contributions from Im(/) to Re(/<(ek,k) = Ret/i 



> and hence of sub-subleading order. This completes the 
(ck, k) -I- (?!)(k) are in total beyond subleading order. 
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F. Re<^(ek,k) to subleading order 

In the following we recover the real part of the gap equation to subleading order by Hilbcrt transforming the 
imaginary part of the gap equation (|31|l and adding the equation for the local gap, (/)(k), Eq. H88() . This shows how 
(f) (j) and 7i[y^pQj|,] ~ cj) combine to a subleading-order contribution. The gap equation for Re0(ek,k) reads to 
subleading order 



Re 0(ek, k) 



3i2n)^ I d^'^^^^Ce..^) tanh ^ g 



E 

(7=± 







^ + ^HDL(ek-0-gq,p) 



?-CI 

J dppTTl{k,p, q) Ag_22(ek - creq,p) 



Tr*(fc,p, (7)A*jDL(ek-CTeq,p) 



(90) 



where we used p [lo, p) = for p > Agi in the effective theory, cf. Eq. IjBSap . Furthermore, all terms ^ coth have been 
neglected. Adding Eq. H88() . the 1/p^-term from the soft electric gluon propagator in Eq. (|90() restricts the p— integral 
of (j) from Agi to 2/i. This is the aforementioned cancellation of (j) and 7i[y^pQjj,], which reduces these terms to the 
order g4>. After approximating the hard magnetic gluon propagator as Aq 22(ek — o'eq,p) = l/p^ + 0(Aq/Agi), one 



can combine it with the remaining contribution from 
arrives at Eq. (124) of Ref. 



Using Tr\{k,p, q) - Tr* (fc,p, g) = 4 + 0(Aq/Agi), one finally 



G. Im (/>(ek + 27?, k) exponentially close to the Fermi surface 

In Sec, nil B1 and ITTTni the contributions A and B to Imt/) have been estimated for different regimes of lo and cf. 
Tab. IIIIVIIII In the case that uo ^ Ky with g < y < 1 and ^ < Ay/^ we found Al^^^ to be the dominant contribution 
to Irntj). The cut of the longitudinal gluons is suppressed by a factor (0/M)*', while the gluon poles do not contribute 
at all. In other regions of lo and C different gluon sectors are shown to be dominant. Furthermore, for lu > 2mg also 
the contributions from B would have to be considered, since there B is suppressed relative to A only by one power of 
g and therefore contributes at subleading order to Im0. For a; ~ Ay with g < y < 1 and C < Ay/3 we find for the 
imaginary part of the gap 

r^e r 2M^ lu* p^ 



Im0(c^ + zry,k) ~ j ^ ^^(eq) Re </.(eqq) J dp ^ 



Ai A 

y 



2 



0/dy'sinf^] fl '"^y 



9(27r)2 \M J ^ J " \ 2 J \ 2 
1 

^ .g0|cos(^)+0(g2), (91) 

where we substituted uj = Ay and d£,/^ — dy' \n{(f>/M) and used ln(0/M) = —3TT^/{^/2g). Furthermore, it was 
sufficient to approximate Tr*(fc,p, g) ~ —2. This result agrees with Eq. (81) in Ref. [13 where a different approach is 
used. 



IV. CONCLUSIONS AND OUTLOOK 



In this work we studied how the non-local nature of the gluonic interaction between quarks at high densities affects 
the energy and momentum dependence of the (2SC) color-superconducting gap function at weak coupling and zero 
temperature. For this purpose, energy and momentum have been treated as independent variables in the gap equation. 
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By analytically continuing from imaginary to real energies and appropriately choosing the contour of the integral over 
energies, we split the gap equation into two coupled equations: one for Ret/) and one for Im0. 

In order to solve these equations self-consistently, the gap had to be estimated for all energies and for all momenta 
satisfying j/c — /i| < Aq, where Aq ~ is the quark cutoff of the effective theory employed in this work. For quarks 
exponentially close to the Fermi surface, we have proven the previous conjecture that, to subleading order, one has 
(j) — Re (j), where Re 4> is the known subleading order solution for the real part of 0, which neglects all contributions 
arising from the non-analyticities of (f). 

Furthermore, we found that, exponentially close to the Fermi surface and for small energies, only the cut of the 
magnetic gluon propagator contributes to Im cj). Thus, the analytic solution of the imaginary part of the gap equation 
is rather simple, cf. Eq. ||MJ. For energies of order rUg, we showed that also the electric cut and the gluon poles 
contribute to Im</), cf. Tables imvml The increase of the imaginary part with increasing energies can be interpreted 
as the opening of decay channels for the quasiquark excitations. The peak Im cj) ~ g^/i occurring for energies just 
above rUg reflects the decay due to the emission of on-shell electric gluons. 

Treating energy and momentum independently, the solution also includes Cooper pairs further away from the Fermi 
surface, up to \k — fi\ ^ gfi. This becomes important when one is interested in extrapolating down to more realistic 
quark chemical potentials where the coupling between the quarks becomes stronger: With increasing g also quarks 
away from the Fermi surface participate in Cooper pairing. These are not included in the Eliashberg theory, where 
one assumes that Cooper pairing happens exclusively at the Fermi surface. 

Finally, it would be interesting to generalize our analysis to non-zero temperatures. The dependence of (j3{T)/(f){T — 
0) on T/Tc, where Tc is the critical temperature for the onset of color superconductivity, agrees with that of a weakly 
coupled BCS superconductor if one neglects Imcp 0|- For strongly coupled superconductivity in metals it is known 
jl9j |. however, that Im cj) is significantly modified at non-zero temperatures due to the presence of thermally excited 
quasiparticles. This in turn gives rise to important deviations from a BCS-like behavior of 0(T) /0(T = 0) at energies 
larger than the gap. An analogous analysis for color superconductivity would be an interesting topic for future studies. 
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APPENDIX A: SPECTRAL REPRESENTATION OF THE GAP FUNCTION 

The function (l){K) which solves the gap equation for imaginary energies ko = i{2n+ l)7rT, must be analytically 
continued towards the axis of real frequencies kg ^ lu + ii] prior to any physical analysis. It is shown in the following 
that the gap function must exhibit non-analyticities on the axis of real energies, in order to be a non-trivial function 
of energy. We require that the gap function converges for infinite energies, i.e., (piK) (^(k) for |fco| ^ oo. Then 
<j>(K) can be written as 

</>(if) = 4>{K) + m , (Ai) 

where the energy dependence of (l){K) is contained in (j>{K) and (f>{K) for |fco| oo. The non-analyticities of 
(j){K) are contained in its spectral density 

P4,{uj, k) = [(l){uj + iTj, k) - (j){uj - iTj, k)] . (A2) 
Zm 

With Cauchy's theorem one obtains for any /cq off the real axis 

oo 

0(fco,k)^ / d^M^ + 0(k), (A3) 
J UJ - ko 

— oo 

where the first term is identified as (j){K) in Eq. l|Aip . For ko = oj + ie one obtains 

oo 

(j>{oj + ie,k) = V I dw' P^Y''^'^ -f0(k)+i7rp^(co,k) . (A4) 

— oo 
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Hence, for real ^^(k) and ^^(wjk) 



Re</>(w + ze,k) = V I dcj' ^^^^^^ + 0(k) , (A5) 



Im + le, k) — 7rp0(cj,k). (A6) 

Furthermore, one finds the dispersion relations for + 177, k) 

00 

Re ^ w + ze, k = -V / dw' — ^^-^ — , A7a 

TT J UO' — UO 

— 00 

oc 

Im + ie,k) = r / div' -i- ^ , (A7b) 

n J LO ~ uj 

—00 

i.e., Im4>{ijj + ie, k) and Re(/)([J + ie,k) are Hilbert transforms of each other, Re^ = 7i[Im(^]. It follows that (f> 
only if both Re (j) and Im (jj are nonzero. Consequently, the gap function (j) is energy-dependent only if it has a nonzero 
imaginary part, Im0 = Imcf) 7^ 0, which in turn is generated by its non-analyticities along the real kg axis, cf. Eqs. 
(|X2|) and l|X6|) . 

The energy dependence of the gap function (p{K) is constrained by the symmetry properties of the gap matrix 
^^(K). It follows from the antisymmetry of the quark fields that the gap matrix must fulfill 

C^+iK)C-^ = [<i>+{-K)f , (A8) 

cf. Eq. (B4) in 22]. Since C^^A^C^^ — ["/bA^y.]'^ and in the 2SC case [J3T2]^ = J3T2 it follows for the gap function 

^{K) = 4>{~K) . (A9) 

Assuming that the gap function is symmetric under reflection of 3-momentum k, (^(fco,k) ~ (^(fco,— k) one obtains 
with Eqs. HA5IA6|I 

Re0(w + i77,k) = Re(f){-uj + iri,k) , (AlOa) 
Im 0(0; + 177, k) — — Im (/)(— + i?/, k) , (AlOb) 
p^{uj,k) = -p^{-uj,k) . (AlOc) 

Consequently, Re (j) is an even function of uj while Im cf) and are odd. In the case of the color-flavor-locked (CFL) 
phase the gap matrix has the structure ^ J • / in (fundamental) color and flavor space, where the matrices 
{Jf)gh = —i^fgh and (j'^)™" = _^^kmn -j^ color and flavor space, respectively 0. Since (J • = J ■ I, Eq. 
(|A10alA10cll are vahd for the CFL phase, too. 

Inserting Eq. (|A10c() into Eq. (|A3(I one flnds that (l){K) is real on the axis of imaginary energies, whereas (/'(-ft') is 
complex for energies off the imaginary axis. 

APPENDIX B: SPECTRAL REPRESENTATION OF THE GLUON PROPAGATOR 

In pure Coulomb gauge, the gluon propagator ||SJ) has the form 

A°"(P) = A\P) , A°^(P) = , A'^{P) ^ {S'^ ~ffP) A*(P) , (Bl) 

where A^'* are the propagators for longitudinal and transverse gluon degrees of freedom. In the spectral representation 
one has [23| 

00 oc 

A^(P)^-1 + / d^^'("'P) , A*(F)= / dc.^*(^'P). (B2) 
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For hard gluons with momenta p > Agi one finds 



while for the soft, HDL-resummed gluons with p < Agi one has j2^ l85l Is^ 

PiA^^P) = Pif^'i^^P) {6 [u; - ivijip)] +6[LU + uJiAp)]} + Pef{^,p)0{p- |w|) , 



where 



UJ {uj^ — p^) 

p2 (jp. _|_ _ ^2^ ' 

2M2 



pT^^'i^.p) = 
pT\^.p) - 



TT p 

UJ (lu^ — p^) 
imi up- — (up- — p2)2 



p + w 



p — UJ 



M2 W p2 



TT p p^ ^ LO^ 



2M- 



UJ 



2 3 2 / ^ 



In 



p + 



p — UJ 



n 2 



(B3a) 
(B3b) 



(B4) 

(B5a) 

(B5b) 
(B5c) 

(B5d) 



APPENDIX C: (NON-) GENERATING THE BCS LOGARITHM 

The mixed scale Aj^ introduced in Eq. (|36|l can be used to analyze the generation of the BCS logarithm. For this 
purpose we have split the integral over ^ = k — fi according to 



M 

2 / 



Al 

2 / 



As 

2 / d^ 



Ao 

2 / dC 



(CI) 



Al 



where we abbreviated 4>q = Re0(eq,k). The first term is readily shown to be of order g'^(f). The last term is of order 
g^(j), which is found after noting that <j)q ^ g cj), cf. Eq. 1351) . and that A^ is smaller but of the order of M. The second 
term can be written as 



9 — 4>q 

J fq 

Al 



g 

^"(^) j'^y sin [^)^gcp 



(C2) 



where use was made of Eqs. (|34I36|) and the BCS logarithm ln{M/(f)) being of order 1/g. It is shown that the BCS 
logarithm arises from integrating over intermediate scales, Ai < ^ < ■'^g- This observation is useful for estimating 
numerous integrals. Note, however, that in the full QCD gap equation one also has the gluon propagator under the 
integral. Then the region Ai < ^ < Ag is additionally enhanced. 

The contributions arising from P4,{uj) in the gap equation are suppressed, because the oddness of p^[u}), cf. Eq. 
(|A10c|) . prevents the generation of the BCS logarithm. At many places, cf. Eqs. (|59I6UI62I65I66I71I77|) . this oddness 
gives rise to logarithmic dependences of the following form 



d£, , 



$ + A, 



Al 



A^ 



(C3) 



where < y < 1. To show that the BCS logarithm is prevented we introduce the dilogarithm [37 





Li2(a;) 



ln(l - , 



(C4) 
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which has the values — j^tt^ = Li2 (—1) < Li2 (a;) < (1) = ^tt^ for — 1 < x < 1. We write the first term on the 
r.h.s. of Eq. (jU^ as 



Ai Ai/A„ 

Since < Ai/A,y < 1 this term is of order one and no BCS logarithm has been generated in this term. The second 
term on the r.h.s. of Eq. I|C3I) is 



Ao Ao/Ay A„/Ao 



Ay 1 1 

In the second step one substituted x = with dx/x — "d^/S,. Similarly to Eq. IjCSp . < Ay/Ao < 1. Hence, also 
this term is of order one and no BCS logarithm has been generated here, either. 
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